Generalized parton distributions in the deuteron 
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Abstract. We introduce generalized quark and gluon distri- 
butions in the deuteron, which can be measured in exclusive 
processes like deeply virtual Compton scattering and meson 
electroproduction. We discuss the basic properties of these 
distributions, and point out how they probe the interplay of 
nucleon and parton degrees of freedom in the deuteron wave 
function. 



Introduction. The partonic structure of the deuteron has 
been explored in terms of the parton distributions ac- 
cessible in deep inelastic scattering [||, and in terms of 
the form factors measured in clastic lepton-deuteron pro- 
cesses [||J|] . It is natural to ask what can be learned from 
generalized parton distributions (GPDs), introduced not 
long ago in . For the nucleon case it has been shown 
that these quantities contain unique information about 
the dynamics of quarks and gluons in QCD bound states, 
beyond what can be unraveled from ordinary parton dis- 
tributions and form factors. Here we extend these studies 
to the case of the deuteron, with the aim of providing the 
theoretical framework to analyze and interpret present 
and future measurements with deuteron targets. We re- 
strict ourselves to parton distributions of twist two, and 
to the parton helicity conserving sector, which is relevant 
in most phenomenological applications. Quark and gluon 
helicity flip GPDs can be treated with the same methods. 

Generalized quark distributions for the deuteron. As in 
the nucleon case, the GPDs for the deuteron are de- 
fined through non-diagonal matrix elements of quark- 
antiquark operators on the light cone. Their general de- 
composition can be written as 
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where n is a light-like four- vector, n 2 — 0. The incoming 
and outgoing deuterons respectively have momenta p, p', 
helicities A, A', and polarization vectors e = e(p, A) and 
e' = e'(p',\'). We write P = p + p' and A = p' - 
p, and choose Ji's variables x, £ = — (A.n)/(P.n), t = 
A 2 as arguments of the GPDs Hi and Hi. The tensors 
and depend on the four-vectors p, p', and n. 
One need only keep tensors which do not vanish when 
contracted with e Q , e'p, given the orthogonality conditions 
e.p = e'.p' = of the polarization vectors. With the 
constraints from parity invariance we find that the 
are linear combinations of the five tensor structures 



{gp a , Ppn a , npp' a , p p p' a , n n a }. 
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Similarly, the A^ a are linear combinations of the seven 
tensors 



{fp/ia/p'", Wa^V' C^oflV", 

P t "p" J n p n a , e^pp p p p' u n p p' a , 
f-nvpa p p p' v n p np, e^ pa p^p lu n p pp}. 
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Using the Schouten identities || one can show that only 
four out of these seven are linearly independent. The 
first x moments of generalized parton distributions are 
elastic form factors. To keep the corresponding relations 
simple we take among the tensors V^ % \ Ay> those which 
appear in the conventional form factor decomposition of 
the vector and axial currents : 

(p',A'|^(0) 7 ^(0)b,A) = -G x {t){^.e)P^ 

+ G 2 (t) [e»(e'*.P) + e'*»{e.P)] 

-G 3 (t) ( e .P)( e '*.P)— , 

(p', X'\m f 75 V>(0) \p, A) = - iGi(t) # afh e'^P^ 

a e 7 fe'* P) + e'* 7 fe P) 
+ iG 2 (t) e p aPl A a pP li! L^f [ > , (5) 

where our convention for the antisymmetric tensor is 
e oi23 = 1 and M is the deuteron mass. The matrix ele- 
ments are here defined flavor by flavor; to get the conven- 
tional form factors, one must weight with electromagnetic 
or weak charges and sum over flavors. For the matrix el- 
ements of the non-local operators we define 
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Since it is determined by the quark operators in Eqs. (|lj) 
and (||), the Q 2 evolution of the GPDs Hi and i?j is 
exactly the same as for spin 1/2 targets, worked out in 

BUS- 

Time reversal properties and sum rules. The action of the 
time reversal operator on the matrix element (|l|) leads to 
the relation 



e >*P e * V (2(pA,n)H(x,tt) 
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where we have made explicit the dependence of the ten- 
sors on the relevant four-vectors. Taking the com- 
plex conjugate of Eq. (Ill) we get 



^e* a vff(P,A,n)H?(x,Z,t) 
= e* V" K« (P, - A, n)JJ, (*,-£,«). 
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For the axial vector case, we obtain relations analogous 
to (0) and (|) by replacing V {l) with A® and H with 
Combining these conditions we find that all nine GPDs 
are real. However, their behavior under time reversal is 
not uniform and we have: 

H i (x,Z ) t)=H i (x,-S ) t) (i = 1,2,3,5), 

Hi(x,S,t)=-Hi(x,-S,t), 

Hi(x,£,t) =Hi(x,-£,t) (i = 1,2,4), 

H 3 (x,t,t) = -H 3 (x,-Z,t)- (9) 

Note that in the non-forward case, time reversal invari- 
ance fixes the phase of the generalized parton distribu- 
tions and determines their behavior under sign change 
of the skewedness parameter £, but does not limit the 
number of GPDs fig ]. Integrating V\>\ and Ay\ over x 
one obtains the local matrix elements (^) contracted with 
//(P.n). Since the tensors that accompany the distri- 
butions Hi, Hi in Eq. (0) are linearly independent, we 
obtain the sum rules 



dxH i (x ) ^t) = G i (t) (i= 1,2,3), 



dxH i {x,£,t) = G i {t) (i = 1,2), 
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dxH 3 (x,£,t) = 0, 



dxHi(x,£,t) = 0. 
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The integrals over P4, P3 and P5, P4 do not correspond 
to form factors of the local vector or axial currents and 
therefore vanish. In the case of H4 and H 3 this is due to 
time reversal constraints, whereas the definitions of P5 
and H4 involve the tensor n^n 1 ' / (P.n) 2 , whose analog 
cannot appear in the decomposition of the local currents 
due to Lorentz invariance. 

The forward limit. Let now study the forward limit of the 
GPDs, which defines the usual parton distributions. In 
the parton model, i.e., at leading twist and leading order 
in a s there are three independent structure functions in 
deep inelastic scattering, Pi, 61, gi, whose probabilistic 
interpretation in terms of quark densities reads Q 
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Here q^^(x) represents the probability to find a quark 
with momentum fraction x and positive (negative) helic- 
ity in a deuteron target of helicity A. The unpolarized 
quark densities q x are defined as q x (x) — q x (x) + q x (x). 
From parity one has q x = q^ X - The densities for anti- 
quarks are defined in analogy. Note that the probabilistic 
interpretation for Pi and gi is similar to the one in the 
spin 1/2 case, whereas the function 61 does not appear 
for spin 1/2 targets. In the forward limit the only struc- 
tures in Eq. (^) that survive are those proportional to 
Hi , P5 and Pi , because in that limit we have A = and 
e.P = e'.P = 0. Using the results for helicity amplitudes 
given below, one gets: 



Pi (x, 0,0) 



q 1 (x)+q- 1 {x)+q°(x) 



H 5 (x,0,0) = q°(x) 
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for x > 0. The corresponding relations for x < involve 
the antiquark distributions at —x, with an overall minus 
sign in the expressions for Pi and H 5 . With Eq. @ we 
thus have 
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and recover the parton model sum rule bi(x) — of 
PH , which was obtained under the assumption that the 
quark sea q — q does not contribute to this integral. 

Helicity amplitudes. In the region £ < x < 1 our GPDs 
can be represented in terms of amplitudes for the scat- 
tering of a quark on a deuteron |10|] , defined as 



Ax>±,X± = - {Vy\ ± A x >\) 



(14) 



with ± referring to the helicities of the quarks. With the 
constraints 
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from parity invariance there are nine independent quark 
helicity conserving amplitudes. Since Ao+,a+ = -Ao-,o- 
we have only four quark helicity dependent distributions 
Hi, compared with the five quark helicity independent 
Hi. Time reversal invariance gives 



A(x, t)xn,x>» = (-1) V X A(x, t)x>»M 



(16) 



and thus does not further reduce the number of GPDs, 
as remarked above. To define the polarization of the 
incoming deuteron we introduce 
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and e(±l) = = F(e ( ' 1 ^ ± ie^)/y/2 then correspond to def- 
inite light-cone helicity [jl2|. This approximately coin- 
cides with usual helicity in frames where the deuteron 
moves fast, provided that sgn(p 3 ) = — sgn(n 3 ). The po- 
larizations for the outgoing deuteron are obtained from 
Eq. (jl7]) by the exchange p^ «-> £ <-> — £ and an 
overall sign change for and e^ 2 ). With this we get 
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where D = (to — t)/(AM 2 ). The remaining amplitudes 
can be easily obtained from the relations (|l^) , (|l^) and 



Ax' —^.x—pi Ax' [i,Xfi(Hi 



■Hi). 



(19) 



As required by angular momentum conservation, one 
gets a factor y/to — t for each unit of helicity flip. 
Note that #2,4 and H4 only appear with longitudinal 
deuteron polarization, and that the only GPDs appear- 
ing in double helicity flip amplitudes are H 3 and H 23 . 
In the forward limit we have A(x, 0, 0)a+,a+ = <7j (x), 
A(x, 0, 0)a-,a- = qf (%) and find the relations (|l2|). 

Gluon distributions. Let us turn to the gluon distribu- 
tions in the deuteron. Instead of the matrix elements ([!]), 
(B) we now have 
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with F a/3 = \e al3 ~* s F lS . We take the same tensors 
yM, as f or q uar k distributions, given in Eq. (||). 
Note that the Hf are even and the Hf odd in x. 
Their behavior under time reversal is the same as 
in (^|) for the corresponding quark distributions, and 
the definitions and expressions of the helicity ampli- 
tudes A 9 ^, A u are also analogous to the quark case. 
The forward limit is now A 9 (x, 0, 0)a+.a+ = xgb[x) and 
A 9 (x, 0, 0)a-,a- = xg^(x), with an extra factor x com- 
pared to the quark case. 

Some Phenomenology. The deuteron GPDs can be ac- 
cessed in hard exclusive processes such as deeply virtual 
Compton scattering in ed — > edj, and electroproduction 
ed — > edM of a meson or a meson pair Jjjjj . The rele- 
vant kinematical limit is that of large invariant momen- 
tum transfer Q 2 to the lepton at fixed xb and t, where 
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the Bjorken variable xb is defined as in deep inelastic 
scattering. Factorization formulae are the same as for 
nucleon targets (|JD1, with the appropriate replacement 
of the hadronic matrix elements, but the same hard scat- 
tering kernels. The Q 2 behavior of the amplitudes and 
the selection rules for photon and meson helicities also re- 
main the same, since they depend on the hard-scattering 
process, not on the target spin. Electroproduction of a 
pseudoscalar meson selects the GPDs Hi, vector meson 
production involves the H i and Hf , and all distributions 
appear in Compton scattering. Notice that the isosinglet 
nature of the deuteron simplifies the flavor structure of 
the GPDs and thus of the scattering amplitudes. One 
consequence is that pion exchange does not contribute 
to any of these processes, in contrast to the nucleon case, 
where it may give important contributions through the 
quark distribution E |14|| . In kinematics where the Bethe- 
Heitler process dominates in electroproduction ed — > ed'j, 
one can use the methods of to study Compton scat- 
tering through the interference of the two processes. At 
sufficiently large Q 2 this interference term gives access to 
a linear combination of GPDs, weighted with the electro- 
magnetic deuteron form factors G\^.z{t)- 

Discussion. To get a feeling for the physics of GPDs in 
the deuteron, consider the approximation where they are 
written as a convolution of the nucleon GPDs with the 
light-cone wave function ijjp+ n for a proton and a neu- 
tron in the deuteron p7| . The struck nucleon then has 
to absorb the entire momentum transfer, in particular its 
plus-component parameterized by £. Assuming for sim- 
plicity that ipp+n is only nonzero if the plus-momentum 
fraction of the proton in the deuteron is between — w) 
and |(1 + w), one finds that all Hi and Hi vanish for 
£ > w. The £ dependence of the deuteron GPDs thus re- 
flects the width of the wave function ip p + n in longitudinal 
momentum fraction. For £ well above w they provide ac- 
cess to deuteron wave function components that cannot 
be described in terms of individual nucleons. To unravel 
such components in inclusive deep inelastic scattering at 
xb > 1 has turned out to be difficult. Approaching this 
region from below, the struck quark has to take the en- 
tire momentum of a nucleon in the convolution picture, 
so that not only the deuteron wave function ip p + n be- 
comes small but also the parton density in the nucleon. 
This is not the case for GPDs with their independent mo- 
mentum variables x and £ : at £ > w one can still have 
any value for the plus-momentum fraction of the struck 
parton in the target. 

We also note that, since in a convolution model b\ 
requires a d wave component in ip p + n j the same holds 
for H5 . Double hclicity flip amplitudes also need a d wave 
admixture: the helicity flip of a nucleon cannot exceed 
one unit, so that orbital angular momentum is necessary. 
Hence H$ and Z/2,3 involve the d wave part of tp p + n - 

In this letter we have focused on the case where the 
deuteron scatters elastically. One may extend our study 
to the case where it dissociates into a proton and neutron, 



or a more complicated hadronic system, introducing ap- 
propriate transition GPDs. Factorization still holds in 
the processes discussed above, provided that the invari- 
ant mass of the dissociative system is small compared 
with the hard scale Q 2 . 

Let us finally stress that for low enough £ the general- 
ized parton distributions for elastic deuteron transitions 
are by no means small. Neither is there any suppressing 
factor in the cross sections. It should thus be possible 
to perform exclusive electroproduction experiments on a 
deuteron target where those on a nucleon are possible. 
We hope for a rich harvest of physics on this topic in the 
forthcoming years at existing facilities such as DESY and 
Jefferson Lab. 
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